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We compute the QCD running coupling on the lattice as defined from the 3-gluon vertex. We present 
the results of an exploratory study at /3 = 6.0 on a 16^ lattice, which show that for momenta larger than 
2 GeV, the coupling runs according to the 2-loop asymptotic formula, allowing a precise determination 
of the A parameter. Our renormalization procedure corresponds to a momentum subtraction scheme in 
the continuum and, most remarkably, one does not need lattice perturbation theory to match the results 

to WS scheme. We obtain a'^'' (Mz) = 0.115 ± 0.003 ± 0.008. 
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We report here the preliminary results of a recent lattice 
determination of the QCD running coupling, a^, from 
the 3-gluon vertex ^. This method differs from the pre- 
vious ones as the coupling is determined in terms of the 
fundamental fields and vertex functions of the contin- 
uum Lagrangian. This is achieved by evaluating two 
and three-point off-shell gluon Green functions, in the 
Landau gauge, and imposing non-perturbative renorma- 
lization conditions directly on them ^. By varying the 
renormalisation scale, the method allows the determina- 
tion of as{q) for different momenta from a single simula- 
tion. Of course this is valid provided that the renormal- 
isation scales can be chosen in a range of values between 
the non-perturbative region and the high momentum re- 
gion, where discretization effects become important. An 
analogous non-perturbative approach has been recently 
applied to the renormalisation of composite fermion op- 
erators ^. As will be explained in the following, a crucial 
feature of this procedure, which corresponds to a momen- 
tum subtraction renormalisation in continuum QCD, is 
that renormalised Green's functions do not depend on 
the lattice regularisation. As a consequence, lattice per- 
turbation theory (LPTH) is not needed in order to relate 

MS 

the measured coupling to , and the relation between 
the two schemes can be entirely computed in continuum 
perturbation theory. Despite recent proposals ■* to im- 
prove the convergence of LPTH series, this step still pro- 
vides one important source of systematic errors in the 

MS 

final prediction for , expecially at low /3 values. Fi- 
nally, given the conceptual simplicity of the method, the 



extension to unquenched QCD presents in principle no 
extra problems. 

2 The Method 

Provided systematic errors are under control, in order to 
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compute from the lattice it is sufficient to: 

1) fix the scale of momenta by determining the lattice 
spacing a; 

2) choose a suitable renormalisation scheme which defines 
the renormalised coupling; 

3) match the result to the MS scheme. 

As for step 1) we take the value of determined by 
Bali and Schilling ® in their string tension measurements, 
that is a-i = 1.9 ± 0.1 GeV at /3 = 6.0. We quote a 
systematic error on the scale to take into account the 
uncertainty resulting from other possible choices. 

In our renormalisation scheme, step 2), the running 
coupling at scale /j, is defined from the renormalised am- 
putated 3-gluon vertex, as 



rg^(p^,o,-p///) 



On the lattice the amputated bare 3-gluon vertex is de- 
fined by evaluating first the complete 3-point function 

(3) 

Gg\^p^{api,ap2,ap2,) and then amputating its external 
legs by dividing it by the gluon self-energies correspond- 
ing to the external momenta. We work in the Landau 
gauge and choose the asymmetric kinematical point de- 
fined by 
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This is a choice suitable for the lattice symmetries, and 
such that the tensor structure of the vertex reproduces 
the tree-level one ^. Then in the a — limit, the bare 
3-gluon vertex can be written as 
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{pa,0,-pa) = 2 Zy'^{pa) goia) pp. 
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where go is the bare coupling constant. Notice that on 

(3) 

the lattice the calculation for Tg' yields directly the 
product Zy^go in a non-perturbative way. The renor- 
malised 3-gluon vertex is than related to the bare one 
through (the a — limit is understood in the following) 
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where Za is the gluon wave function renormalisation. 
In order to compute Za , we impose that for a fixed mo- 
mentum scale p^ = fi"^ the renormalised propagator takes 



its continuum tree-level value, 
condition is simply 



Thus the renormalisation 
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where (pa) is the scalar part of the bare gluon self- 
energy in the Landau gauge. The above equation defines 
Za in a non-perturbative way. 

In this way one can express gRijJ-) entirely in terms 
of two and three-point bare Green's functions non- 
perturbatively evaluated on the lattice 



, , \pf G^^lp„ipa,0,-pa) ^ 

9R(I^ — 7, ; \p2 = n^ 

2pp ^(2)^^^. f^) 



(6) 



The above formula is the basis of our computational pro- 
cedure. 

3 Numerical Results 

The lattice calculation was performed on a 16K Con- 
nection Machine 200. We generated two data sets, each 
comprising 150 gauge configurations, on a IG'* lattice, at 
j3 = 6.0 and 6.2. The [3 = 6.2 lattices were produced 
essentially for testing purposes, so that in this paper we 
only quote numerical results obtained at /3 = 6.0 A cru- 
cial step in the method is the accurate implementation 
of the lattice Landau gauge condition by using a Fourier- 
accelerated algorithm. 

The calculation of the renormalised running coupling 
then proceeds as follows. As an initial step, we evaluate 
the gluon propagator and compute Za from eq.(5). Next, 
we compute the complete three-point function G g of the 
gluon field. Finally, g{lJ.) is obtained from eq.(6). This is 
plotted as a function oi jj, = \p\ in Fig.(l). 
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Figure 1: Running coupling vs. momentum. 

The tree-level form of a bosonic lattice propagator 
suggests that higher order terms in the lattice spacing are 
small when the lattice momenta satisfy p = 2Tr n/ La <^ 
T^I'la, where n is an integer. As our lattice size is _L = 
16, we limit lattice momenta values to n < 3 in each 
direction. In addition, we will also discard momenta such 
that p'^ > (Gn/La)^. The important physical question 
is whether there exists a range of momenta for which 
the coupling that we measure runs according to the 2- 
loop asymptotic formula. To answer this question, we 
compute A.~--, as a function of the measured values of 
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g(l^) according to 
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where bo = ll/167r2, bi = 102/(167r2)2 and A^ 

the QCD scale parameter for the renormalization scheme 

that we are using. In the asymptotic region, we expect 

A.^r^-:, defined from the above equation to be a constant. 
MOM ^ 

We plot ^j^f^j]^ versus q in Fig. (2); it appears that 
for q > 2 GeV the data are consistent with a constant 
value for the QCD scale parameter. In order to extract 
from the above data a prediction for A.~--,, we fit the 

MOM] 

data points for /i > 2 GeV to a constant. This yields 
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0.96 ±0.02 ±0.09 GeV, 



(8) 
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where the first error is statistical and the second one 
comes from the uncertainty on the value of a~^ . 

4 Matching to WS 

In order to extract from the numerical results for A 
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a prediction in the MS scheme - step 3 in section 2 
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estimate (Mz). One can take into account, at least 
partially, the effect of quenching, by assuming that the 
quenched and the 3-flavor couplings would be approx- 
imately the same at scales as low as ~ 0.7 GeV. A 
different procedure to estimate the effect of quenching 
is described in ref.^. By averaging the two results and 
evolving the coupling up to the Mz scale, we would ob- 
tain 



{Mz 



0.115 ± 0.003 ±0.008 



(12) 



where the first error is dominated by the uncertainty on 
, and the second one keeps into account the residual 
effect of the quenched approximation. This result is com- 
patible with the first unquenched calculations, see ®, and 
with other quenched lattice computations ^. 

5 Conclusion 



Figure 2: A parameter vs. momentum. 



only need the ratio ^Jfs/ ^mom' which can be computed 
entirely in continuum perturbation theory. 

The general relation between the A's in two schemes 
A and B, can be written as 
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+ 0(g'M). (9) 



In our case, in the Landau gauge and in the quenched 
approximation, = 0, we get 
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0.35 



(10) 



Note that the scheme that we have adopted in thejion- 
perturbative calculation differs from the usual MOM 
scheme as it contains an extra constant term in the ver- 
tex renormalization constant ^ . This is a linear term in 
the momenta, not proportional to the tree level, equal on 
the lattice and in the continuum, which is perturbatively 
included in (10). Using the nonperturbative result (8), 
we get 
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340 ± 10 ±40 MeV 
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This is the main result of our computation. We empha- 
sise that even for /3 as low as 6.0, the uncertainty from 
the one-loop matching condition is smaller than the er- 
ror on the determination of the scale and comparable to 
the size of the statistical error. Of course this depends on 
the fact that we only use continuum perturbation theory, 
and it would not be true with other methods which rely 
on LPTH. 

At this point we face the problem of extracting a 
prediction for the full theory. In particular, we aim to 



We have shown that a non-perturbative determination 
of the QCD running coupling can be obtained from first 
principles by a lattice study of the triple gluon vertex. 
Systematic lattice effects appear to be under control in 
our calculation. In ref.^ a comparison with a correspond- 
ing LPTH calculation provides some important consis- 
tency checks. In the future we plan to study further the 
role of infrared and ultraviolet cutoffs, by increasing the 
physical volume and the range of /3 values for our lattices. 
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